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The theory of two counter-propagating polarized beams interacting in a 
randomly birefringent fiber via the Kerr and Raman effects is developed and 
applied to the quantitative description of Raman polarizers in the undepleted 
regime. Here Raman polarizers, first reported by Martinelli et. al. [Opt. 
Express. 17, 947 (2009)], are understood as Raman amplifiers operating in 
the regime in which an initially weak unpolarized beam is converted into an 
amplified fully polarized beam towards the fiber output. Three parameters 
are selected for the characterization of a Raman polarizer: the degree of 
polarization of the outcoming beam, its state of polarization, and its gain. All 
of these parameters represent quantities that are averaged over all random 
polarization states of the initially unpolarized signal beam. The presented 
theory is computer friendly and applicable to virtually all practically relevant 
situations, including the case of co-propagating beams, and in particular to 
the undepleted as well as the depleted regimes of the Raman polarizer. © 
2010 Optical Society of America 
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1. Introduction 

Fiber optic Raman amplifiers constitute an integral part of many contemporary high-speed 
optical networks. Owing to their broad amplification bandwidth, Raman amplifiers suc- 
cessfully compete with erbium-doped-fiber amplifiers, p. By the nature of the process of 
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stimulated Raman scattering (SRS) in silica fibers, Raman amplifiers are highly polarization- 
dependent: the gain coefficient reaches its maximum value whenever the state of polarization 
(SOP) of the signal is parallel to the pump SOP. Whereas the amplifier gain is reduced by 
two orders of magnitude if the two SOPs are orthogonal, [2H1]- In real fiber networks the 
fiber birefringence changes stochastically along a fiber span. The net effect of this stochasti- 
cally varying birefringence is to greatly reduce the polarization dependence of Raman gain. 
As a result, the effective Raman gain reduces to its span-averaged value, which is equal to 
approximately half of its maximum value. Such averaging of the gain between two orthogonal 
states of polarization is a favorable feature, since polarization-dependent effects are undesir- 
able in most of the telecom applications which have been developed so far. The value of the 
polarization-mode dispersion (PMD) coefficient D p quantifies the degree of gain averaging, 
which appears to be rather efficient for D p > 0.2 ps/y/km. For smaller values of D p , the 
onset of polarization-dependent Raman gain is typically combated by scrambling the pump 
SOP. 

In this paper our interest diverges from the common trend, as it focuses on the exploita- 
tion of the polarization dependence of Raman gain. Motivated by the recent developments of 
transmission protocols based on polarization multiplexing [5] , we are interested in developing 
nonlinear devices for achieving all-optical and ultrafast polarization control. Thus, we are 
interested in the use of fiber with low values of the PMD coefficient, and do not perform the 
scrambling of the pump beam. Moreover, the stronger the polarization-dependent gain, the 
better the performance of our device - a Raman polarizer. These devices are Raman ampli- 
fiers that are being fed by a weak unpolarized light signal, and convert it into a powerful yet 
highly polarized light towards the output. Preliminary theoretical investigations of Raman 
polarizers and their proof-of-principle experiment, both in the co-propagating geometry, have 
been recently reported in Ref. [6]. This geometry has been analyzed in more details in the 
undepleted regime in Ref. [7] 

Here we set our goal in developing a simple and physically transparent theory of Raman 
polarizers which is applicable to both co-propagating and counter-propagating geometries, 
and suitable for further use as a tool for analyzing the operation of these devices. Our goal is 
also to address the two main issues which are related to the performance of Raman polarizers. 
The first issue is the degree of polarization (DOP) of the outcoming signal. Given that the 
average DOP of the scrambled signal is zero initially, we say that our goal is reached if the 
DOP of the output beam is close to unity. In this case, we shall refer to such Raman amplifiers 
as ideal Raman polarizers. The second issue is provided by the output signal SOP; namely, 
its relation to the pump SOP. The DOP and SOP of the outcoming signal fully characterize 
the performance of Raman polarizers. The third important quantity is the value of gain of 
the Raman polarizer. We shall comment on this issue, too. 



2 



It is important to distinguish from the outset between conventional Raman amplifiers 
and Raman polarizers. Typical Raman amplifiers for telecom applications are long, often 
surpassing 10 km, require pump powers of the order of 1 W, and operate with relatively 
high values of the PMD coefficient. In contrast, Raman polarizers are shorter, requiring a 
fiber length of the order of 2 km, pumped with high power light - of the order of 5 W and 
higher, and favor low values of the PMD coefficient. The last two characteristics substantially 
strengthen the polarization dependence of the gain, as we shall see from the subsequent 
analysis. Therefore the operational regimes of Raman amplifiers and Raman polarizers are 
drastically different, thus preventing the results of the theory of Raman polarizers to become 
a trivial copy of the results already known from the theory of Raman amplifiers. 

The description of polarization-dependent gain in randomly birefringent fibers is based on 
the vector theory of SRS, as it was developed for instance in Refs. [8HTT]. Theories presented 
in Refs. 0,[TO] bring analytical insights into the problem, but their applicability is limited 
to the regime in which the beat length of the fiber Lb is substantially smaller than the 
birefringence correlation length L c : Lb <C L c . Because of this restriction, those theories are 
not suitable for the proper description of Raman polarizers, for which the opposite inequality 
usually holds, see Ref. [7] and the analysis that follows. In contrast, the theory in Ref. [IT] 
is most general and accurate, however it is resource-consuming in terms of computational 
time, as it requires 10 3 ~ 10 4 longer integration time than the theory which is presented here. 
Moreover, since the theory of Ref. [H] is not formulated in terms of deterministic differential 
equations, it is rather difficult to extract from it physically transparent information about 
the role of the different processes ruling the overall polarization dynamics. 

As a matter of fact our theory is essentially the generalization of one-beam theory of Wai 
and Menyuk in Ref. p2] to the case of two beams interacting not only via the Kerr but also via 
the Raman effect. It is close to the approach that was undertaken in Ref. [10], but it proceeds 
till the end with virtually no approximations. The only important requirement here is that 
min(L, Lnl) 3> L c , where L is the total length of the fiber and L^l its nonlinear length. 
This inequality holds true for almost all practically relevant situations. Another assumption 
in our theory is that in the randomly varying birefringence tensor characterizing the fiber 

AB = A/3(w) (cos 6a 3 + sin 6tr x ) , (1) 

we shall only treat the orientation of the birefringence axis 9 as a stochastic variable, while 
we keep the value of the birefringence A/3 as a constant (this is the so-called fixed modulus 
model). It is now a well established fact that the theory with 9 as the only stochastic variable 
and the theory where both 9 and A/3 are stochastic variables (random modulus model) both 
produce virtually identical results. In Eq. ([[]), <j\ and <x 3 are usual Pauli matrices. Note that 
we define the birefringence beat length as Lb = 27r/A/3. 
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2. Model 

Starting from the most basic model of interaction of two fields within a Kerr- and Raman- 
active medium, whose tensorial response as it is relevant to fused silica is properly taken into 
account, see Ref. [9], we follow the lines of derivations that were outlined in Ref. [13] and 
formulate the equation of motion of the pump beam in the form 

c)X J c^X/J 1 

±i-^r + WM-fiT = ~ abu p - p PP W; ■ u p )u p + (u p ■ u p )u;} 

-\i PS w: ■ u s )u p + {u s ■ u p )u: + (u p ■ u:)u s ] + ie p9o (u p ■ w s )u s . (2) 

Here U p = (U px , U py ) T is the field vector written in terms of normal polarization modes e x 
and e y . f3(u) p ) is the propagation constant and (5'{oj p ) is its frequency derivative. The upper 
sign ("+") describes the configuration when the signal and pump beams propagate in the 
fiber in one direction (co-propagating geometry), while the lower sign ("— ") is related to the 
situation when they propagate in opposite directions (counter-propagating geometry). The 
self-polarization modulation (SPM) coefficient 7 PP is the usual nonlinear coefficient of silica, 
7, calculated at the central frequency u p of the pump beam. Cross-polarization modulation 
(XPM) coefficient ^ ps is also equal to j(co p )] g is the Raman gain coefficient; e p = —u p /u s . 
The equations for the signal beam and all coefficients are the same as above but with indices 
p and s interchanged, and with e s — 1. 

First, let us transform the equation (|2]) to the local axes of birefringence by performing 
the rotation of the field vector Uj as 

Oi=f cosf 9 sin !W (3) 

V ~ Sm 2 COS 2 / 

where j = p, s. Equation (T2]) is not altered by this transformation, the only difference being 
the change of the form of the birefringence tensor, which now becomes 

Here, 9 Z is the derivative of 9 with respect to z. It is different from zero owing to the random 
changes of orientation of the birefringence axes. Namely, the change of 9 is driven by the 
white noise process: 9 Z = ge(z), with zero mean and (g e (z)gg(z')) = 2L~ 1 S(z — z'), where L c 
is the birefringence correlation length, as mentioned in the Introduction. 

In line with our assumption that min(L, Lnl) ^> L c , we separate the fast motion related to 
the rapid changes of 9 from the slow motion that occurs on the scale of the nonlinear length. It 
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is important to note that this separation is exact and yet does not involve any approximation. 
The approximation is to be made at a later stage. The required transformations read as 

V p = T p (z)U p = ( a \ a l)u p , (6) 

Here V p = (V p i, V p2 ) T and V s = (V s x, V S 2) T ■ Matrices T p and T s obey equations of motion 

QT 

±^+A5( Wp )T p = 0, (8) 
QT 

i—^ + AB(cu s )T s = 0, (9) 

oz 

These transformations eliminate the birefringence terms from the equations of motion of V p 
and V s and bring about a vast number of cubic terms composed of different combinations of 
VpU Vp2, V s i, V S 2 and their complex conjugates. Factors in front of these terms are products of 
two coefficients of the form u m u n , or u* m u n , or u* 2 u* n , where m, n — 1, . . . , 14. Products with 
m = n we shall call self-products, while with m ^ n cross-products. Here, ni = |ai| 2 — l a 2| 2 , 

a\bi — a 2 &2) u 8 — —(pia 2 + &2 a i); u 9 = ^(^i a 2 — 0*^2)' n io = — i{a*bi + 02^2)5 n n = a i^2 + bia 2 , 
U12 = aib\ — a^2? M i3 = —i(aibi + 0^2)5 u u = «(ai&2 — °2^i)- 

In the thus obtained equations of motion for V p and V s we perform the ensemble average 
(over realizations of the random process that describes the linear birefringence). Thus, we 
write (u m u n ) instead of u m u n . This change holds true only in the limit when the stochastic 
variations are faster than the nonlinear beam evolution. This is exactly the place in the 
derivation where our single approximation comes into play. At this point we also need to 
apply the ergodic theorem 

{/)= lim- f 'dz'f(z'). (10) 
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Our goal is to calculate ensemble averages of all necessary self- and cross-products, and on 
this way we may complete the derivation of the differential equations for V p and V s . 

The equations of motion for u n with n — 1, . . . , 14 can be easily formulated basing our- 
selves on equations ([8]) and (Q. As these equations are linear, in order to find an ensemble 
average of any function of these coefficients ( in our case pair products) we need to construct 
a generator. We refer to the Appendix in Ref. [12] for details of this procedure, and only give 
here the final result. With this generator we are able to formulate the equations of motion for 
the ensemble averages of the products of the coefficients. Thus the solutions to the equations 
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of motion 



5* = O, (id 

BG 2 2 



Bz L c 
dG 3 



(G 1 -G 2 ) T 4A/3K)G 4 , (12) 
±4A/3(a; p )G 4 , (13) 



«/~< 

_i = __ G4 ± 2A/3(^ P )(G 2 - G 3 ) (14) 

yield the result for the self-products {(uf), (^ 2 )}, 
{(Re 2 ( M4 )), (Re 2 (M 5 )), (Re 2 K))}, 

and {(Im 2 (ii 4 )), (Im 2 (ti 5 )), (Im 2 (ti 6 ))}, if we associate them with {G±, G 2 , G 3 ] with initial 
conditions given as (1, 0, 0), (0, 1, 0), and (0, 0, 1), respectively. 

The remaining self-products {(Re 2 (u 7 )), (Re 2 (tig)), (Re 2 (ug), (Re 2 (uio))}, 
{(Im 2 (u 7 )>, (Im 2 (u 8 )), (Im 2 (u 9 ), (Im 2 (u 10 ))}, 
{(Re 2 (u n )), (Re 2 (u 12 )>, (Re 2 (u 13 ), (Re 2 (u 14 ))}, 

and {(Im 2 (wn)), (Im 2 (u 12 )}, (Im 2 (u 13 ), (Im 2 (u 14 ))}, can be found from the equations 

BC 2 

— (Gi - G 2 ) + 2 [±A/?K) - A/?(u; s )] G 5 , (15) 



Bz L 
dG 2 _ 2 

dz ~ L 
0G 3 

dz 
BG 4 



- r (G l - G 2 ) - 2 [±A/3(w p ) + A/%,)] G 6 , (16) 

2[±A/3(u; p ) + A/3(u; s )]G 6 , (17) 

-2[±Af3(u p )-Af3(u s )}G 5 , (18) 

— ^ = [±A/3(u; p ) - A/3(w a )] (G 4 - d) - — G 5 , (19) 
oz L c 

= [±A/3(u p ) + A/3(u s )} (G 2 - G 3 ) - -G 6 , (20) 

when we associate them with {G±, G 2 , G 3 , G 4 }, with initial conditions as (1, 0, 0, 0), 
(0, 0, 0, 1), (0, 1, 0, 0), and (0, 0, 1, 0), respectively 

In order to find the cross-products we constructed appropriate generators and found that 
all the cross-products that are of interest to us turn out to be equal to zero. Similarly, terms 
of the form Re(w n )Im(w n ) also vanish. Thus, many SPM, XPM, and Raman terms in the 
final equations of motion disappear. The thus found equations of motion for the fields are 
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conveniently formulated in Stokes space. They read as 

U^+ftuJ^jS® = 1 (u p )^xjf(z)S^ + S^xJ x (z)S^ 

+e p g [S { S) J R0 S^ + S®J R (z)SV) , (21) 

(JL + ^ = l{Us) ( S(S) X J S ] ^ S{S) + sM x J x(z)S {p) ) 

+9o (sIpJmSM + S { s) J R (z)S^) . (22) 

Here is a three-component vector: S^ = (S^, S 2 \ S3), whose components are defined 

as 

S( = VjiV* 2 + V; x V j2 , (23) 

si = KVflVp - v 3l v; 2 ) , (24) 

^ = l^i| 2 -|^ 2 | 2 , (25) 

where j = p, s. It is instructive to formulate separate equations for the powers of the signal 
and pump beams, which may be defined through the components of the Stokes vectors as 

S^=(s? 2 + Sf + Sf) 1/2 : 

e P 9o (JroS^S^ + J R1 (z)S[ s) S[ p) + J m (z)sPs® + J m (z)Sj?>S®) , (26) 




9o (JroS^S^ + J R1 (z)S{ s) S[ p) + J m (z)SPsP + J m (z)S ( 3 s) S^) , (27) 

Matrices in equations (I2ip and fl22|) are all diagonal with elements 
J R = diag(J m , J R2 , Jjh), J x = diag(J X i, J X 2, Jxa), Js = diag(J S i, J S 2, Jsa)- These ele- 
ments are different for the counter-propagating and the co-propagating interaction geome- 
tries. 

In order to complete our theory, we need to express all elements in these matrices in terms 
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of ensemble averages of self-products: 



Jri — 


(Re(w 2 4 - u 2 10 )) , 


(28) 


JR2 = 


-(Re(u 2 u + u 2 w )) , 


(29) 


JR3 = 


/ 1 1 2 i 1 2\ 
-(|«14| - \Ul0\ ) , 


(30) 


,J A 1 


2 

— ('Refu?r, -1- u 2 n — u 2 , — it, 2 *)) 

XJ-^V^IO ~ 13 9 "14// ' 


(31) 


Jx2 


2 

— CRe(ll 2 r> + 111 A — II 2 . — 1l 2 n)) 

<j \ uc rlO T "14 u 9 13// ' 




Jx3 = 


9 

/I |2 , 1 |2 1 |2 1 |2\ 
g\|(X9| -r | t*i4| \ a l3\ \ a 10\ ]i 




Jsi = 




(34) 


Js2 = 


-^(Re(^)>, 


(35) 


JS3 = 


I [3<t^> - 1] , 


(36) 



and also Jr = (|^io| 2 + |^i4| 2 )- 
3. Physical considerations 

Model equations (I2ip and ( )22|) represent the central finding of our paper. These equations 
describe vector Raman amplification under rather general conditions. The SPM, XPM, and 
Raman matrices contain z-dependent elements on their diagonals, whose values are expressed 
in terms of the self-product coefficients in Eqs. fT28|) -f l36|) . while the coefficients themselves 
are found as the solutions of equations (1TT ]) -(!H ]) and (fT5"|) - (l2"Ul) . 

We are about to apply these model equations to the problem of Raman polarizers. Before 
treating this specific case, some comment about the Raman matrix is in order. It is commonly 
accepted that the net effect of polarization-dependent gain is the attraction of the signal SOP 
to the pump SOP. While this is certainly the case for isotropic fibers, an extension of this 
idea to the case of randomly birefringent fibers should be undertaken with caution. As a 
matter of fact, we found that in general in randomly birefringent fibers the signal SOP is 
not attracted to the pump SOP. In order to demonstrate this, the most straightforward way 
is to analyze equation f )27|) for the power of the signal beam. For isotropic fibers, we get 
Jri = Jr2 = Jr3 = 1, so that parallel signal and pump SOPs provide the maximally possible 
gain. 

When we turn to the case of randomly birefringent fibers, we find that the elements of 
the Raman matrix are not always equal to each other in the co-propagating geometry, as it 
was demonstrated in Ref. [7]. Moreover, these elements are never equal to each other in the 
counter-propagating geometry, as it is exemplified in Fig. [U As can be seen in this figure, the 
Raman matrix elements are not only different from each other but they can also be negative. 
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Fig. 1. Elements of the Raman matrix Jr = diag(J#i, J^ 2 , Jm) (Jri - black 
solid, Jr2 ~ red dashed, and J# 3 - green dotted) as function of distance in the 
fiber for Lb{oj p ) = 50 m and L c = 1 m. (note that the black solid and green 
dotted curves nearly coincide, i.e. Jri = Jr%. 



Although in general the signal SOP is not attracted to the pump SOP, the analysis of the 
co-propagating configuration in Ref. [7] shows that in the regime when Raman amplifiers act 
as ideal Raman polarizers the output signal SOP is always almost perfectly aligned with the 
pump SOP. For this regime, the polarization properties of the randomly birefringent fibers 
are virtually the same as in the case of isotropic fibers. 

On the other hand in the counter-propagating geometry the situation is totally different. 
In the regime when Raman amplifiers act as ideal Raman polarizers, the typical relation 
between the diagonal elements of the Raman matrix is as shown in Fig. [TJ Thus, the first 
and the third components of the signal Stokes vector are attracted to the corresponding 
components of the pump Stokes vector, while the second components are repelled from each 
other. Therefore whenever the pump SOP contains an appreciable admixture of the circular 
polarization, the signal SOP is never attracted to the pump SOP. In the next section we 
suggest a simple rule with a wide range of applicability on how to determine the output 
signal SOP. 

The polarization-dependent gain is quantified by the values of the diagonal elements of the 
Raman matrix. The larger these values, the better the performance of the Raman polarizers. 
In fibers with high PMD values the diagonal elements quickly vanish near the input fiber 
end already. For smaller PMD values the ^-dependent elements of the Raman matrix keep 
appreciable values across the whole fiber span, as demonstrated in Fig. [U Indeed, low PMD 
fibers are good candidates for the implementation of Raman polarizers. 
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The elements of the XPM matrix exhibit a similar dependence upon the PMD coefficient 
and have a magnitude which is comparable to that of the diagonal elements of the Raman 
matrix. Nevertheless the nonlinear polarization rotation which is due to the XPM interaction 
is very weak, and it has virtually no effect on the performance of Raman polarizers operating 
in the undepleted regime. 

The SPM effect for the signal beam has also no impact on the performance of Raman 
polarizers. First, this is because the diagonal elements of the SPM matrix Jg^ vanish on 
first 100 m of the fiber for the practically relevant range of beat lengths and correlation 
lengths, a range that we define as follows: 



which is given here in km. Second, because the signal beam is too weak to experience a 
significant nonlinear self-interaction, especially near its input end. 

In contrast, the SPM effect can in principle be sizable for the pump beam. Given that the 
pump power is relatively high (8 W in our simulations), even the first 100 m are enough to 
perturb the pump SOP. However, these perturbations remain relatively small (of the order 
of 0.1%) for a pump power as high as 8 W. 

4. Numerical results 

We performed extensive simulations in a configuration where the pump and signal beams 
counter-propagate through a 2 km long fiber, and the pump power is set to 8 W at z = L. The 
input power of the signal beam is set to be equal to 0.1 mW, which is low enough to ensure 
that Raman amplification occurs in the undepleted regime for all z and the entire range of 
fiber parameters. We varied the beat length and the correlation length within the range that 
is defined in Eq. f[3"T|) . We followed two quantities of interest at the output end of the fiber at 
z = L as function of the beat and the correlation lengths. The first such quantity is the DOP 
(D) of the signal beam, which is the central characteristic of Raman polarizers. Starting 
with an unpolarized signal beam (with D = 0), we aim at finding such parameter regimes 
for which the outcoming signal beam has DOP close to unity. The unpolarized signal beam is 
modeled in the Stokes space as an ensemble of beams with N ~ 5000 states of polarizations 
that are uniformly distributed over the entire Poincare sphere. On the other hand the pump 
is a polarized beam whose SOP is varied as detailed below. The second quantity of interest 
is the so-called alignment parameter, which is defined as 



where the angle brackets indicate average over all N realizations of the SOPs of the input 
scrambled signal beam. This quantity characterizes the average orientation of the signal SOP 



0.001 < Lb < 0.1 and 0.0001 < L c < 0.05 , 



(37) 




q{s) q 

°0 °( 



(38) 
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Fig. 2. DOP of the signal beam (black, solid) and alignment parameter Af± 
(red, dotted) as function of correlation length L c for four pump SOPs: a) 
3- 1 /2(i ; i ) i). b) (1, 0, 0); c) (0, 1, 0); d) (0, 0, 1). The value of the birefrin- 
gence length Lb(ou p ) is indicated on the plots in km. Two circles on plot d) 
indicate one (of infinitely many) pair of points with equal PMD coefficients. 
Other parameters are: input signal power 0.1 mW, input pump power 8 W, 
Raman gain go = 0.6 (W-km)" 1 , nonlinearity parameter 7 = 1 (W-km) -1 , 
linear fiber losses a = 0.2 dB/km, total length of the fiber L = 2 km. 



with respect to the pump SOP. Note that in the copropagating geometry the appropriate 
alignment parameter reads as [7], 

(s[ s) s{ p) + s { 2 s) si p) + s ( 3 s) si p) ) 

A n = ± rv -^ '- . (39) 

I I q {s) nip) v ' 

This alignment parameter is nothing but an average cosine of the angle between the pump 
and signal SOPs in Stokes space. 

We start with the characterization of the performance of the Raman polarizer in the 
counter-propagating geometry by following the D and for four different SOPs of the 
pump beam and two values of the beat length, as a function of the correlation length. The 
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Fig. 3. DOP of the signal beam for two pump SOPs as function of the correla- 
tion length L c for the value of the birefringence length a) Lb = 0.016 km and 
b) Lb = 0.05 km. Other parameters are the same as in Fig. [2j These two pump 
SOPs are those which maximize (black, solid) and minimize (red, dashed) the 
DOP of the output signal beam. 

corresponding results are shown in the four panels of Fig. [2j All four plots demonstrate that 
small correlation lengths and large beat lengths favor the efficiency of Raman polarizers. In 
the first three plots the maximum DOP is 0.94, while on the fourth plot it is of only 0.84. 
Such large difference in the maximal achievable DOPs gives a hint that the performance of 
Raman polarizers is rather sensitive to the choice of the pump SOP. 

How sensitive is this dependence is demonstrated by the two panels of Fig.[3j For generating 
each point on these plots we took 256 realizations of the pump SOPs, and for each of this 
SOP we followed the signal DOP as a function of the correlation length. We plot only those 
two realizations which provide maximum and minimum DOP. From the analysis of these 
plots we may conclude that from the viewpoint of efficient Raman polarizers only a limited 
range of pump SOPs leads to a good performance, while the rest of pump SOPs perform 
poorly. The particular SOP which maximizes the DOP depends on both the correlation and 
the beat length values. 

In all of the four plots of Fig. [2] we may observe that the alignment parameter is relatively 
close to unity whenever the DOP is more than 0.9. From this result we draw the conclusion 
that on average the SOP of the outcoming signal beam is related to the pump SOP as dictated 
by the alignment parameter that was introduced in Eq. fl38|) . A more detailed information 
on the signal SOP is available when basing ourselves on the following definitions of the 
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Fig. 4. Alignment factors defined in Eqs. (HDJ) - (J42J) : A$ (black solid); A$ 
(red dashed); (green dotted), as function of the correlation length for 
L B = 0.05 km and for two values of the pump power: a) 6 W; b) 8 W. 
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which are calculated at the Raman polarizer output z = L. These expressions quantify the 
pairwise proximity of the components of the signal and pump Stokes vectors. The dependence 
of the alignment factors upon the correlation length for Lb = 50 m and for two values of the 
pump power is plotted in Fig. HI One may observe that the output signal SOP depends on 
the pump power. Therefore, signal polarization stabilization may be deteriorated by relative 
intensity noise of the pump laser source. 

Another important issue is the selection of fibers for Raman polarizers. Usually, this selec- 
tion relies on the value of the PMD coefficient. In most practical situations, this coefficient 
fully characterizes a randomly birefringent fiber. However, in some cases the knowledge of 
this parameter only is not sufficient for making a conclusion on the performance of Raman 
polarizers, a fact that was first noticed in Ref. [7]. For example the plot in Fig. EJ^d) demon- 
strates that two fibers with equal PMD coefficients lead to Raman polarizers with rather 
different performances: in one case D m 0.26, while in the other case it is 0.1 only. In the 
regime of ideal Raman polarizers (with DOP> 0.9) and in the considered parameter range, 
the PMD coefficient is indeed a reliable characteristic, in the sense that fibers with equal 
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Fig. 5. Gain of a Raman polarizer as function of the correlation length. The 
state of polarization of the pump beam is (1, 0, 0), the signal beam is unpo- 
larized. Values of the beat length Lb are indicated on the plot in km. Other 
parameters are the same as in Fig. |5J 

PMD coefficients provide Raman polarizers with similar performances. Nevertheless, it is 
often desirable to consider the beat length and the correlation length separately, as we do in 
this study, rather than unite them under the single PMD coefficient. Note that for L ^> L c , 
which is always the case for Raman polarizers, the PMD coefficient is defined as in Ref . [12] , 



Finally, let us consider the gain characteristics of Raman polarizers. Owing to their strong 
polarization-dependent gain, the gain of a Raman polarizer is sizably larger than that of 
typical polarization insensitive Raman amplifiers (i.e., with high PMD values or that use 
polarization scrambled pump beams), see Fig. |5j This increase in gain demonstrates that 
Raman polarizers are simultaneously very efficient Raman amplifiers. This feature means 
that the amplification of signal beams with powers in the mW range will almost certainly 
lead to the depletion of the pump. This regime is left aside in this study, as our goal was the 
demonstration of the basic properties of Raman polarizers only. The regime with depleted 
pump is so rich and parameter-sensitive that it becomes rather difficult to draw general 
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conclusions. However, if necessary our theory is capable of treating the depleted regime too. 

5. Conclusions 

We have developed a general theory for describing the interaction of two beams in randomly 
birefringent fibers via the Kerr and Raman effects. This theory can be applied to both the 
co-propagation and counter-propagation configurations. Here we considered only the counter- 
propagating geometry, while details on the co-propagating regime can be found in Ref. [7]. 
A comparison between these two studies shows that the counter-propagating case is more 
demanding, in that it requires for instance PMD coefficients below 0.008 ps/ ^/km for getting 
DOP slightly above 0.9 in a 2 km fiber. Whereas the co-propagating geometry can provide 
DOPs as high as 0.99 for a D p as large as 0.014 ps/ Vkm with only a 1.5 km long fiber. In 
both cases the pump power is 8 W. 

However, it is important to point out that the counter-propagating configuration has the 
advantage of clamping the signal SOP to the SOP of the pump source. In the co-propagating 
geometry the output signal SOP is also related to the pump SOP, but both depend on the 
stochastically changing birefringence of the fiber, thus complicating the control over the 
output signal SOP. 

In addition to the theoretical model of the Raman amplification in randomly birefringent 
fibers, we presented the scheme for the quantification of the performances of Raman polar- 
izers. We have identified three main characteristics of Raman polarizers: the DOP of the 
outcoming signal beam, its SOP defined in relation with the pump SOP, and the amplifier 
gain. The present study was limited to the undepleted regime only, although the theory is 
readily applicable to the depleted regime as well. 
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